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GRAVITY  turbulence  connected 
WITH  INTERLACES 


C.  M Tchen 
The  City  Col iege  of 

The  City  University  or  New  York,  New  York 
ABSTRACT 

The  spectral  distributions  of  turbulence,  as  generated  by 
the  gravity  waves  on  the  interface  between  two  fluids,  are  investi- 
gated. Both  stable  and  unstable  su' faces  are  considered.  An  un- 
stable surface  refers  to  the  early  development  of  turbulence  from 
the  Taylor  instability  A stable  surface  may  refer  to  the  sea 
surface  A repeated  cascade  method  is  used  to  close  the  hierarchy 
of  correlations  at  their  fourth  order,  and  to  describe  the  eddy 
transport  property  as  a memory  chain  or  eddy  relaxations.  The 
production,  inertia  and  dissipation  subranges  of  spectral  distribu- 
tions on  an  unstable  surface  with  friction  are  found  to  follow  the 

-2-3  -3  -35-1 

laws  k"  , k , and  k~  for  the  k.net.c  energy,  and  k , k , and 

_5 

k for  the  surface  elevation  The  inertia,  eddy  dissipation  by 
gravity,  and  molecular  dissipation  subranges  on  a stable  surface 

_3 

with  friction  are  found  to  follow  the  law  k for  the  kinetic 

-1  -5  -5 

energy,  and  the  laws  k , k and  k for  the  surface  elevation, 
respectively  The  spectra  with  surface  tension  are  also  investigated. 
The  physical  parameters  and  the  numerical  coefficients  are 
determined  analytically 

-I- 


For  neutral  fluids,  we  can  distinguish  hydrodynamic  turbulence 
and  gravity  turbulence  fhe  i at  ter  moy  occur  in  a free  medium, 
eg,  in  atmospheres  and  oceans  with  a mean  gradient  in  density  or 
temperature  It  can  also  appea'  on  a stable,  or  an  unstable, 
surface,  of  a heavy  1'quid  bplow,  or  above,  a light  one.  For 
example,  tne  sea  surface  's  such  a stable  surface  The  turbulent 
motions  in  a staole  background  may  require  an  external  source  of 
energy  for  the' r ms'ntenance  On  the  other  hand,  an  unstable  sur- 
face may  eventually  develop  into  singular  ringers,  broken  bound- 
aries and  droplets  which  .cease  to  constitute  a continuous  surface. 
Therefore,  the  gravity  turbulence  can  only  refer,  to  the  early  stage 
of  the  deve'opment  of  a continuous  turbulent  surface  under  an 
unstable  cond it  on 

Mathemat  >ca  I ly  spe.-ixing,  when  the  pressure  is  eliminated 
between  the  Nav ie' -Stokes  equation  and  the  equation  of  continuity, 
the  hydrodynam  ^ turbulence  s de:- . - bed  by  one  single  equation  of 
motion,  and  the  gravity  turbulence  described  by  a system  of  two 
equations,  governing  the  velocity  and  temperature,  or  density,  in  a 
free  medium,  or  the  velocity  and  the  elevation  on  a moving  surface. 

In  view  of  the  complexity  of  the  dynamical  equations,  the 

dimensional  method  had  been  rei  ed  upon  for  solutions.  It  enabled 
I 2 

Kolmogorort  and  He  senberg  to  derwe  the  spectral  laws  for  the 
velocity  fluctuat  ons. 


r(L  ) - - t i A)  - ■ 


(i) 
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in  the  inert’a  and  diss  pation  sub^-anyes  of  a hydrodynamic  turbu- 
lence On  the  same  bas's,  Shu  ^ prcpo’-ed  the  spectrum 

i~  - co^„t  /V^  £ J i 


tor-  the  gravity  turbulence  n a t-ee  medium,  and  Phillips 
obtained  an  energy  spectrum 


/ = 4 "k  ~ 

tf 


and  the  elevation  spectra 


ti  (/&)=■  'i. 


;-~j  fu)j  - ^0 


in  gravity  turbulence  with  :urfac?s  The  spectra  may  be  determined 
in  the  space  of  wave  number  k,  or  o*  frequency  w The  parameters 
in  the  above  torwuiai  are:  the  rate  of  energy  dissipation  e, 

the  kinematic  viscosity  , tne  acceleration  of  gravity  g,  and 
the  Brunt- Vaisa ia  frequency  N representative  of  a mean  gradient 
in  temperature  or  density 

it  is  obvious  that  the  spectra  should  not  be  identical 
fo''  both  stab  e and  unstable  surra  .es  Since  the  dimensional 
analysis  on'jr  recogn.zes  the  pa  ameters,  without  going  into 
dynamical  mechan i sms , it  is  not  able  to  determine  the  conditions, 
stable  O'  unstable,  under  wh^ch  the  dimensional  laws  (2),  (3)  and 
(4)  should  apply,  although  the  authors  intended  for  their  use 
under  a stable  condition 

In  view  of  the  above  difficulties,  it  is  necessary  to  con- 
sider an  analytic  treatment  sine  most  analytic  theories  in 
hydrodynam’C  turbulence  cannot  even  predict  satisfactorily  a 


I 


Kolmogoroft  law  (1)  which  requires  only  one  parameter,  we  shall  resort 

5 

to  the  method  of  repeated  cascade,  proposed  by  Tchen  . It  enables 
closing  the  hierarchy  ot  cor  re ! at  ions  to  their  fourth  order,  and 
characterizing  a turbulent  transport  property  by  means  of  a memory 
chain  ot  eddy  relaxations  For  a hydrodynamic  turbulence,  it 
derived  the  Kolmogoroft  law  of  turbulence  in  the  inertia  subrange, 
and  a k 1 law  in  a gradient  flow  In  the  present  paper,  we  shall 
extend  the  method  to  the  gravity  turbulence  on  stable  and  unstable 
surfaces-  The  gravity  turbulence  in  a free  medium  will  not  be 
treated  here 


li  DERIVATION  OF  THE  EQUATIONS 
rOR  THE  MOTION  Or  THE  SURFACE 


The  equat  ons  ot  Nav iC'-Stokes  and  of  continuity  for  the 
motion  or  a liquid  are 


f- 


> a 

! «* 

/ 


f 


(5) 


u. 


(6) 


where  v^t,  Xj , x^,  x^)  is  a veloct>  vector  with  components 
in  three  direction^  x)5  x^,  x^,  p is  the  pressure,  p is  the  density, 
, is  a frict.or.al  coefficient,  and  g is  the  acceleration  of 
gravity 


We  introduce  a velocity  potential  <j>( t , X| , x^,  x^)  such 


that 


[/  « 


“-s  ,M 

i_p_ 

~iX 


/,  Z,  3 


(7) 


and  transform  Eqs  (51  and  (6)  into  a Bernoulli  equation  and  a 
Laplace  equation, 


C 


respect  ve’y,  where  t ( t ) ’s  a function  ot  time,  which  may  be 
determined  Dj  a cond  t>on  at  the  surface,  it  necessary. 

In  order  to  apply  Eqs  (8)  and  (9)  to  the  surface,  we 
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put  X3  -•  ,(t,  Xj , x^)  m (8j,  and  write  the  velocity,  potential 
and  pressure  on  the  surface,  as  fol 'owSi 


h *> ' v)=  v R 


y x * t x 


i SJ 


C#  x, . xj  , l=LZ  , 


f ( r'  x-  XJ  - >z  , / 

R */  ' *z) -- id(£'  *>'  *z  > >i  - £ y - d o-^stctn£  f 

transforming  (8)  into 

(ifj  *■  r<V+  <'  - f 

'J-? 

.Ir-If1;  <-??*•/?  = f(t), 

P f f 

where  T is  a surface  tension 

It  can  be  noted  that  and  f,  which  appear  in  (13), 
can  be  eliminated  by  taking  a gradient,  giving 


I -y t j 

*-  7 -4-  * - 


2 V 

L 


I ^3  (t  \ ' Xi  = *^J 


r vv 

; c,  i 


v 

"2>x  ' ~&y 

t. 


Equation  (14),  combined  with  the  following  kinematic 
equation  for  the  surface  elevation  l 


i_L  _ f'L 
°it  ~ w 


f-  <t  7/  C 


forms  two  fundamental  equations  determining  the  motions  of  the 
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surface  Here 


'id.  -?/<<*, 


(16) 


The  linrts  of  integration  (c,*>)  refer  to  an  unstable  surface,  i.e, 
with  £,  = -1,  and  the  limits  (-*>,;)  refer  to  a stable  surface,  I.e 
with  4 = -1  Equation  (15)  may  be  regarded  as  originated  from  the 
integration  of  an  equation  of  mass  conservation. 

Since  certain  transformations  are  necessary  to  simplify 
Eqs,  (14)  and  (15),  we  write  the  potential  in  fourier  transform 


. -to 


Jjj  (17) 

“UG 

where  the  wave  number  < satisfies  the  relation 


K-2  i-  fc  + 

I A 


- 0 


(18) 


as  a consequence  of  the  equation  of  continuity  (9). 

As  a surface  wave  decreases  its  amplitude  vertically  like 


**/*(-$  i*j ) ' 

we  find 


fc.  » t i (19) 

with  a real  k as  a consequence  of  the  incompressibility  of  fluid, 
as  expressed  by  the  relation  (18).  this  permits  rewriting  (17)  in 
the  form 
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and  reducing  formal  . '/)  to 


rhose  roMTial  -e- at  ions  provides  the  possibility  ot 


tro-. storming  a ■ • guantit'es  re'aied  to  velocity  and  potential 


the  indices  refer  to  1 


Furthermore,  the  functions  R,,  R0,  R,  and  R,  are  corrections 


proportional  to  ; and  its  derivatives,  which  in  turn  are  proportion- 
al to  k*u(k) 

We  ->hai*  now  resort  to  the  BouaSinesq  approximation,  well 
known  in  prob'ems  of  grav i t>  waves , wh  ch  neglects  such  terms  as 
proportional  to  ; or  to  k-u,  except  the  term  gv;  associated  with 

<v. 

the  gravit>  g as  a driving  rorce  .n  the  momentum  equation 
That  approximat  on  amounts  to  neglecting  all  the  above  mentioned 
correction  functions,  and  thereby  reducing  (14)  to  the  followina 
simpler  form  in  x-space  , 

l*T  + ~ v)~  ' ~v (fc)  - H'J  > <29) 

or,  in  Fourier  space. 


^(k-k)«.jk)  <■  rujk) 
~oj~  J 'v  jf  J - ~ v ~ 

L -bo 


~ I 


«*( k}LZ.(k)  = Ejk)  , 


(30 


where 


f(*)  ; (•  - — r 

1 ft 


(31) 


(32) 


and  E is  a driving  force  due  to  the  gravitation  pull  of  the 
surface  The  same  approximation  reduces  (15)  to  the  following 
form  in  Four  er  space' 
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(33) 


Hit! 

~*t 


t>o 


dk 


V c Ay 

4 8 


- £ c-jT  -i  - fAj 

/ 4 ~ 


Here  k - (k^,  k^)  is  a wave  number  vector  in  two  dimensions;  the 
dependence  ot  u and  w’th  t me  is  understood  In  the  following 
we  shall  keep  the  dynamical  eHuat>on5  to'  the  surface  in  their 
Fourier  form,  because  tre  right  hand  side  of  (33)  does  not 
possess  a simple  form  or  inversion 

Here  and  in  the  following  we  omit  the  writing  the 
variable  t in  the  argument  of  ail  time  dependent  functions, 
except  when  a need  tor  spec  fic  distinction  arises. 

Fo'  convenience,  we  define  a speed  of  propagation 


(34) 


and  a dritt  ve'oc'ty  of  the  surface,  Called  "potential  drift," 
m k-space 


w 


A. 


(35) 


With  those  notat’ons  (34;  and  135),  we  transform  ^30)  and  (33) 
to  a more  symmetry  form. 

— — - f y u (k)  +■  { dk  ' I.  k.  Ct.  f k-  k'J  uj  k'J 

oc 

~ iC  ^ -r(U  - ^ (L)  (36) 
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The  form  represented  by  (36)  and  (37)  has  the  advantage  that  the 
gravitational  coupling 


r(i)  ™ E(k)  u.(-k) 

~ E (~k)  u-(k) 

A#  »*  A.  /%» 

is  ot  equal  magnitude  but  opposite  signs  in  the  equation  of 
evolut’-on  ot  kmetic  energy  and  potential  energy 


U,( k)  ■+■{-*■)  y yw(^)  (39) 

A/  -V  ^ A<  /V  ''r 

for  the  case  ot  a stable  configuration,  i.e  4 » -1, 

In  application.',  we  shall  use  both  systems  of  equations, 

Eqs  (30)  and  (33),  or  Eqs  (36)  and  1 37 ) The  letter  system 
will  be  used,  when  the  mechanism  ot  turbulence  is  controlled  by 
the  grav'tat  onal  coupling  between  the  kinetic  energy  and  the 
potential  energy,  while  the  former  system  will  be  used  in  other 
mechan ■ sms  where  such  a coupling  does  not  come  into  play. 
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U:  OUTLINE  Of  THE  SINGLE  AND  REPEATED 


CASCADE  METHODS  Of  TURBULENCE 

1 : 1A  Purposes  of  Single  and  Repeated  Cascade* 

In  the  past,  analytical  theories  of  turbulence  have  been 
based  upon  .nvestigat  ng  the  hierarchy  ot  velocity  correlations 
of  var  ou>  orders,  and  design  ng  the  methods  ot  closure-,  The 
degeneration  of  correlations  from  h’gh  into  lower  order*  can  only 
be  performed  in  an  ambiguous  way  , since  ail  velocities  in  the 
high  o'der  correlation*  have  to  oe  treated  uniformly.  It  is  well 
Known  that  such  a p-actice  may  'ead  to  serious  difficulties 
That  tho*e  velocities  appear  uniformly  is  a certain  consequence 
of  the  mathematical  model  adopted  Physically,  since  a velocity 
consist*  of  many  scale*,  it  ’*  expected  that  like  scales  will 
group  together  to  form  correlates  of  non-un , form  pairs  playing 
non-un’form  physcai  ro'es  In  crde>  to  allow  such  a selection 
ot  scales  in  coo e'ations , we  prescribe  rank*  in  the  velocity 
variable 

A turDulent  motion  * a qua* i -stat ■ onary  process,  having 
a cont  nuou*  spect-um  ot  coupled  scales  The  larqe  scales  form 
a "macroscopic  background,"  prescribing  the  background  conditions 
for  the  mot  on  ot  small  scales  The  smaller  scaie*  move  more 
"randomly,"  and,  as  a result  of  the  statistcai  effect  ot  their 
fluctuations,  shape  up  t'ansport  propert  es  in  the  background 
medium  The  above  division  nto  macroscopic  and  random  variables 
are  relative  to  any  wave  number  o*  the  spectrum  Thus  we  write, 
tor  a veiocty  u m the  phys-cai  space  or  in  the  wave  number  space. 
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{' 


becomes  a "repeated  cascade  " More  yeneraily,  we  can  decompose 
u into 


y. 


(43  j 


instead  ot  (40) , with 


\/W  , 

~ lA.  + +60 

'V  *V  v 

t-)  w W 

v V-  -t  +■  u. 


(44) 

(45) 


1 1 ) ' 

Evidently,  v - u 

/v 

The  superscripts  denote  the  ranks,  with  a higher  rank 
having  a h gher  degree  ot  ' andurnness  Thus  the  ranks  uJ  and 
will  be  cons’dered  as  macroscopic  backgrounds,  while  o' 
and  w i be  random  fluctuations  'n  (40)  and  (43) 

The  rank  u is  »e  ponsible  tor  forming  an  energy,  of 
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scales  up  to  a wave  number  in  tne  spectrum  The  velocity  tr 
represent*,  all  the  ^ m 3 1 i e v scale*  which  shape  up  an  eddy  vis- 
cosity in  the  background  medium  for  the  purpo*e  of  further 
determining  the  app'oach  to  equii.bnum  of  tne  transport  property, 
u’  has  to  oe  subdivided  nto  higher  dnks  v4<)  in  a repeated 
cascade,  where  u^lf  and  u ' ^ ' contribute  to  an  energy  and  a 
re'axation  frequency  required  to  rormrate  the  eddy  viscosity 
In  continuing  the  sequence,  tne  higher  «ank*  u^',  u'^  relay 
to  nigh  cder  relaxation  frequencies  forming  a memory  chain 
The  method  of  closing  the  hierarchies  ar  . sing  from  velocity 
correlation*  of  various  orde»s,  and  from  the  said  memory  chain 
ha*  been  treated  by  Tcheo  'n  applications  to  hydrodynamic  turbu- 
lence  for  the  similar  problem  of  closing  the  hierarchies  and 
the  memory  chain  in  gravity  turbulence,  we  will  need  to  consider 
a repeated  cascade  for  the  variable  u,  while  a single  ca*cade 
suffices  for 


(46) 


and 


W f-  i-v 


(47) 


In  tne  above  p cture  ot  qua* '-stationary  turbulent  proc- 
cess,  the  large  scale  motions  are  considered  relatively  macro- 
scop c and  tne  -ma  ler  scales  a^e  more  'anoom  Thus,  we  can 
associate  a high  degiee  of  randomness  to  high  wave  numbers,  by 
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writing 


yw 


A*),  s 

U.  (xj 


oik  £<,  (k)  J^h  ( i k * ) 

'S  /V  / ^ 


(48) 


-VS 


( ik  xj  > 


(49) 


^ (*J  — j 4 k ix-(k-)  4#  h (c  £,  x ) 

J^-:)  * ' (50) 

-■  I .(;•  ' V k ) -*"j>  (i  k x ) 

J -Co  "/  <v 

In  the  notat  on  (49),  u'J^(k)  understood  to  be  truncated 
between  tne  wave  numbe-  mtervai  (k^-1’1^,  k^)  The  trunca- 
t on  (48)  need;,  not  be  sharp,  and,  t necessary,  can  be  regulated 
by  a scaling  di.-tribut  on 

1 1 1 B Averaging  Rules 

A "cascade  ensemb'e  average,"  o*  "rank  average,"  denoted 
by 

\ > (51) 

is  expected  to  separate  tne  two  components  in  (40),  by  averaginq 
o*er  ■ ea  izations  under  ider.t  cal  mac'oscopiu  background  u 
After  such  an  averaging  procedure,  the  random  component  u1  be- 
comes macroscopical ly  negl  sible,  a^d  the  macroscopic  component 


5- 


u comes  out  intact 

-v 

. / c 

~ U.  , 


Similarly  a 'ank  average 


(53) 


Thus  we  have 

<U  > =*  u,  , <yt  / « O (52) 


would  annui  u , i.e 

a/ 


/ 

v. 


(54) 


and  could  correspond  to  a spatial  average  ot  a length  interval 
X It  irioy  oe  noted  that  x my  tend  to  ; ntimty  in  a homogeneous 
turbulence 

for  the  use  of  distinction  or  a more  general  rank  a 
rank  average 


/ > 

/ 


V ) 


(55) 


corresponding  to  a length  .nterval  X'J^  is  introduced  The 
following  averaging  rules  apply: 


/ 


O 


C/5) 


if 


>s  of 

(56) 

P < yf 

(57) 

as  a generalization  ot  (52)  and  (54). 
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.n  Ref  5 we  have  introduced  a distribution  function  of 


f 

] 


many  velocities  ordered  m ranks  as  random  variables,  and  dis- 
tribution runct  ons  of  reduced  velocity  ranks  They  are  similar 
to  the  distribution  functions  ir.  the  BBuKf  hierarchy  of  kinetic 
theory  or  gas  The^  serve  to  define  the  rank  averaye  (55)  We 
shall  not  ente'  'nto  such  a detail  here 


IliC  RanKs  and  Rank  Values  of  Correlations 

in  v ew  of  the  condition  or  guasi-stationari ty  of  turbulence 
of  rank  cs,  we  can  write 


/.M)  r i , 
\UL  { l>. 


''■A*1;/  X 


(k  i / 

\ v- 


u.^(k')fJd(k^k")  „ (58) 


V 


where 


( */X  7 


(59) 


is  called  <j  “scaling  lactor,"  and  s - 2 in  two  dimensions, 
giving  the  following  re  at'on  between  tne  velocity  correlations 
mi  k-space  and  *-space: 


/ W/  . 

'A. 


/ 

j 


fio 

l I 

<Ak 

->o 


X 


X) 


'(*>) 


X) 


(60) 


In  pd't  Cu la- , we  have 


(6 1 ) 


. 


il 
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reducing  the  t me  Integra'  • 62 ) tu 


It  the  turbulent  mot  on  or  r^n k , ti  of  a 


where 


outfit  enu>  itnd ! 1 scale,  the  aasumpt  on  of  isotropy  can  be  applied 


is  eaMed  dn  edd>  viscosity  ot  the  j-th  rank 


In  view  ot  tne  a^eraong  ruie  (5?)  and  ot  the  expected 


p-eaence  ot 


we  deem  tnat 


in  the  equation  descr'omg  the  evolution  ot  u 


As  a consequence 


as  obta  ned  by  a time  integration,  which 


the  eddy  viscosity 


will  nave  a rank  va iue 


amounts  to  a smooth. ng  process 


lower  Thus 


a '■ank  value  o 


ha>  a rank  value 


A dynatm ca l equaton  tor  the  t’me  evolution  of  a rank 


may  contain  ce  ms  0 


However,  on  account  ot  the  properties 


with  a f 


such  terms  cannot  contribute  to  calculating  statistical  quantities, 

such  as  <uH  u(“M°\  by  the  dynamical  equation.  Thus,  we  have 
J 


< 


u. 


= 0.  if  £<« 


(72) 

(73) 


where  i and  j are  arbitrary  indices.  On  the  other  hand,  if 

6 > ct,  u<*>  raisesthe  rank  value  of  u^  u^  above  a,  so  that 
~ i j 

its  a-th  rank  average  vanishes 

0*X(?O  ^ a y 

\^i  ^ / = 0 / if  ft  • 

For  their  usefulness  of  contributing  to  statistical 
quantities,  we  conclude  that  all  the  terms  in  the  dynamical 
equations  of  single  and  repeated  cascades  must  have  their  rank 
value  conserved.  Nevertheless  a dynamical  equation  for  the 
evolution  of  a rank  a may  still  couple  to  other  ranks,  pro- 
vided they  possess  the  same  rank  value  a. 


For  the 

above  reason,  we 

can  rewrite 

(70)  and  (71)  as 

u it  > - 

\y  v > , 

rank  value 

*-/ 

(75) 

/v  /V  / 

W 

QV  i 

rank  value 

$ <*-2  . 

(76) 

'V 

' V 
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IV.  DYNAMICAL  EQUATIONS  IN 


CASCADE  REPRESENTATION 


As  mentioned  in  Section  III,  we  will  need  a repeated  cas 


cade  for  the  velocity  and  only  a single  cascade  for  the  surface 
Therefore, we  write  the’r  respective  dynamical  equations  as 


V LANGEViN  EQUATION  AND  ONSAGER'S 
RELATION:.  FOR  IURBJlENT  i-lUXE$ 


It  to  be  remarked  that  D,  Dt  represents  a Lagrangian 
derivative,  1 e a rate  ot  change  follow  ng  the  path  of  a fluid 
element  The  variable  t ’n  D Dt  can  be  treated  as  a one-dimen- 
sional  variable  in  the  lay.angian  representaton,  notwithstanding 
its  four  dimensions  m t me  and  three  wavenumbers  in  the  Eulerian 
representat  or.  Under  ;uch  a circumstance,  and  in  analoqy  with 
the  Brown  ar.  movements  ot  molecules,  Eq5  ( 77),  >v  7 8 ) , (82) 


and  (83; 


_an  be  regarded  as  Langevin  equations  for  turbu- 


lent motion,  f k is  taken  n be  a parameter  Tho^e  Langevin 
equations  w'll  be  useful  to  ca'-cuiate  Mu*es  and  aosociated 
transport  coefficients  ot  tu'bulence,  a,  ’t  snould  be  recalled 
that  any  transport  coefficients  are  ndeed  <epresented  by  a 
Lagrangian  descript'on  n statist. cal  thermodynamics  For 
this  purpose,  we  make  a formal  >nteg>ation  ot  (78),  giving 


k)  = flit'  4?\t  £)  -tyj  ■ 


where 


)V*kty*  (’ry 


It  is  to  be  noted  that  a term  s m tar  to  that  in 


1 on  the 


right  nand  side  ot  (82)  ha3  riot  been  car-ied  Over  to  (78)  and 
(85),  for  the  reason  mentoned  in  Il.C 
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Since  a transport  property  contributed  by  a correlation 
from  a rank  v ' in  a qua>  ■ -stationary  background  of  rank  /a^(t,k) 

■v  v v 

in  the  cascade  (43),  the  upper  limit  t will  belong  to  the  quasl- 
stat'onar^  fme  sea  e which  is  much  large'  tnan  the  duration  of 
that  correlat'on  Therefore  that  upper  limit  can  be  replaced  by 

y i \ 

*>,  ana  V 'i.f,k  ) can  be  replaced  by  \r  (t,k' ).  For  the  same 
reason,  the  mtioi  value  wi'i  not  be  cor>elated  with  any 
fluctuat-on  at  t me  t,  thus  simplifying  (84)  to 


t i)  i 4t'  {'  fa',4,)  t-t'fj 


I d V 


(86) 


Here  - t-t 1 , and 


it  \ k")\ f(ilV+E*hl) 


(87) 


We  have  also  neglected  the  f notion  as  being  smail  compared  to 
the  eddy  mining  p'oee^s 

The  enpres-  on  (86)  for  the  fluctuation  avails  itself 

to  formulate  a flux 
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(88) 


All  the  terms  have  the  same  rank  value  s- 1 , except  the  last  term 
on  the  right  nand  sde  of  '85)  which  has  a rank  value  up  to  a-2. 
Such  a dispa  ate  rank  w 1 ■ not  Contribute  to  tlje  flux,  and  there- 
fore wiil  be  omitteo  When  we  make  use  ot  me  property  (58)  and 
the  definition  (54),  we  -•educe  the  riux  to  the  form 


(f't  k-r: 


?-£ 


(89) 


Thus  we  find  the  stat  stma!  etrect  ot  fluctuations  of  rank  u‘ 
upon  the  evolution  ot  o to  take  the  torm  of  a flux  (89),  which 
is  proportion^  to  the  background  velocity  gradient,  with  the 
proportional ity  coeft'c  ent  g-ven  by  the  eddy  viscosity.  That 
a flux  1.'  proport ’ona i to  the  gradient  ot  the  quantity  to  be 
transported  seems  tj  ta*i  unoer  the  general  Onsager  relation  in 
thermodynanr c.-  ot  irrevers.ble  processes  Therefore,  by  repeating 
the  method  used  to'  the  de-ivat  on  of  tne  momentum  flux,  we  obtain 
the  similar  relations  to'  the  su-toce  fluxes 


V.  i ; i f)  - - j , ; k) 

(?.  ( ' '■  ■ t ' > k -k  ;j  . 


(90) 

(91) 


The  relations  1 89 ; , 90;  and  (9i)  wi ! 1 be  called  the  Onsager 

relations  ot  tu-bulence  We  shall  not  nave  the  need  of  fluxes 
of  higher  ranks  a-izing  from  1 ^ orid  , and  therefore  they 

are  not  written  expi t'y  here 
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VI.  TRANSPORTS  Of  MOMENTUM  AND  ENERGY 


By  relying  upon  th  Onsager  relation  (89),  we  can  trans- 
form the  equation  (77)  tor  the  momentum  transport  into  the  form 


where 


\(tk)  = E*t(i,k)  -i U-(t,k)  j 


(92) 


(93) 


is  a relaxation  frequency. 

We  can  assume  an  isotropic  turbulence  for  the  small 

eddies  contributing  n-  , and  write 

J ^ 


(94) 


If  we  treat  (92)  as  a Langevin  equation,  and  integrate, 

we  find 


j ./r  k)  XKjO  (-Cj'v) 


(95) 


Here  the  initial  value  is  omitted,  because  it  does  not  contribute 
in  any  correlation  In  addition  we  have  neglected  the  friction  y 
as  compared  to  the  eddy  relaxation  frequency  u>'. 


— — 


In  a similar  way,  we  reduce  the  equations  (80)  and  (81) 
for  the  transport  of  ;°  and  w°  into  the  following: 


[ 

,jt  + uW] 

tft.k) 

■=  i k ' k u°(t.k ) 

(96) 

*ru; 

* £ LC  k u. °(t  k)  / 

\ v ^ 

(97) 

or,  in  terms  of  E°, 

[it  + a(v 

1 1 

<-v 

= % CZ  It.  k u°(t,k)  . 

(98) 

The  equations  (92),  (96)  and  (97)  for  the  transpors  of 
uc,  ;°  and  wc  take  now  a form  considerably  simpler  than  the 
original  form  in  Eqs  (77),  (80)  and  (81)  The  simpler  form 
casts  its  nonlinear  transfer  into  the  relaxation  frequency  w1. 

Upon  multiplying  (92),  (96)  and  (97)  by  u?(-k),  ;J(-k) 
and  w°(-k),  respectively  and  taking  an  average,  we  find  the 
equations  of  energy  in  Fourier  form: 


Q-tcJ'j  ) u'(-k)}  = i°(l ) E°(-k)  - v u"(k)  U*(- k)  (k+  - k) 


(99) 

(100) 
(101) 


In  the  above  equations  (99),  (100)  and  (101),  the  convection 
terms  do  not  contribute  to  the  energy  evolution  in  homogeneous 
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turbulence,  and  are  therefore  omitted.  The  complex  conjugate 
part  is  represented  oy  (k  - - k),  as  obtained  from  replacing  k by 
-k. 

When  we  multiply  (99),  (100)  and  (101)  by  a scaling  factor 
x°,  and  integrate  with  respect  to  k,  as  prescribed  by  (60),  we 

/v 

obtain  the  equafons  of  energy  balance 


r°-  t*-  */'-  f 


w 


’(-t'f-  f r - t>  <»/; 


with  the  following  transport  functions: 
(a)  for  the  kinetic  energy 


002) 
0 03) 
0 04) 


f-  F(k') 

- iy 

f*  2.Y /V  r(k') 

- Ypk'  X' p(i,k)  «°( i.-kp 

- cq 

r".  fjk  rW  , rt'k)-  y 

005) 
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(b)  for  the  surface 


I 


T 


s 7 { 

rV k'zH(k-)  - Uk'k,iX(sXt,k)t:Tt-k)> 


J 3 2. 


-no 


r;  - j>  r, :W  r(kh 


(106) 


(c)  for  the  potential  energy 

r = y r 

k/  / 


J„=z/ oik’  kiZ Cf(k)  - 1 kiZ  x y(t,-kj) 

*t\  J-  Oo  **  * 


(107) 


In  the  functions  (105),  (106)  and  (107),  we  have  not  written 

out  the  complex  conjugate  part  explicitly,  as  we  did  in  (99), 

(100)  and  (101),  because  the  integration  with  respect  to  dk  in 

the  interval  (-<»,  °°)  rules  out  this  necessity.  The  functions 

T°,  T°  and  T°  ire  called  transfer  functions,  they  govern  the 

cascade  transfer  of  energy  across  each  individual  spectrum.  The 

terms  vJ°,  XJ°  and  xJ°  are  dissipation  functions,  proportional 
c w 

to  the  vorticity  functions  J°,  J"  and  J°,  with  the  molecular 
viscosity  v and  diffusivity  \ as  proportionality  coefficients. 
Although  the  dynamical  equations  (102),  (103)  and  (104)  were 
established  without  consideration  of  dissipations,  the  dissipations 
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are  now  added  to  keep  a conservation  of  energy,  as  an  often 
allowed  practice  in  hydrodynamical  problems.  Finally  r°  and  r° 
are  called  gravitational  exchange  functions,  and  they  govern  the 
exchange  of  energy  between  two  portions  of  separte  spectra  which 
are  subject  to  the  gravitational  pull. 


■ ..«ar r 


VII.  GRAVITATIONAL  EXCHANGE 


The  gravitational  exchange  function  (105)  involves  the 
coupling  <E°(t,k) •u°(t,-k)>° , which  can  be  calculated  by  means 

/V  \>  Ay 

of  (92),  and  gives 


<\Ejt,k)'  u(trk))  = pt  (E  '(t  V £ °(ir-  r>  ~ t1))  ■ (108) 


The  correlation  under  the  integrand  in  (108)  has  its  development 
governed  by  Eq.  (98),  which  is 


(109) 


The  isotropic  form  is  valid  for  the  present  isotropic 
turbulence 


f BZcj,Z  <u/  (i , kj  E°({  -t,  - k))>  , (110) 


and  can  be  integrated  to  give 

(Ea(ttk)  E‘(i.zri)y  - (fftk)-  (-Oi'x) 


, (HD 


where 


0 m kc/z.u)'  <$C 


(112) 
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is  the  ratio  of  the  frequency  of  oscillation  of  a sinusoidal 
gravity  wave  to  the  turbulent  frequency.  Since  the  aeneratinq 
gravity  wave  has  a frequency  lying  in  the  production  subrange, 
and  the  turbulent  frequency  is  in  the  inertia  or  higher  frequency 
subranges,  that  ratio  is  a small  quantity. 

When  we  substitute  the  expression  (111)  for  the  E°- 
correlation  into  (108),  we  find 


r6*  X 


t <Aj  lh  Nt'  E°(t-T  ~k.p 


¥o  Jo 


The  order  of  the  double  integration  can  be  inverted  to  read 


fco  C>o 


dj  cfx  ' jl«I>  f- 2 to ' ■„)  • = J dr  ' J (Jr  ^ fo  (-Zco'z) . 


o r' 


and  reduce  ( 1 1 3)  to 


\EJlV  u = -C(t''~P 

+ 4-  £ d'cj  ' j Jt'  <u k)-£'(t -r  -kj) 


it  is  to  be  remarked  that  the  intearal  term  on  the  rioht 
1 2 

hand  is  of  the  order  of  <ED(t,k)*u°(t,-k)>°,  and  is  a factor 

C.  A,  /V  rs.  /v 

e2  smaller  than  the  term  on  the  left  hand  side.  When  such  a term 
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HHMMI 


is  neglected,  we  simplify  (115)  to 


(116) 


It  follows  the  expressions  for  the  gravitational  exchanges  written 
in  terms  of  the  spectral  function  H(k): 


r 


-do 


(117) 


f . fit ' [t 

_ [k'c(k')f  H(i'J 


u’(k') 


(118) 
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VIII  EDDY  VISCOSITY  AND  MEMORY  CHAIN 


VIIIA.  Structure  of  Eddy  Viscosity 

A tjrbulent  motion,  when  decomposed  into  a macroscopic 
background  and  a random  fluctuation,  finds  an  interplay  between 
the  two  parts  in  such  a way  that  the  statistical  effect  of  the 
random  fluctuations  upon  the  macroscopic  background  appears  as 
a turbulent  stress,  expressed  by  the  Onsager  relation  (89)  as 
derived  by  the  Langevin  equation  (78)  of  fluctuations.  That  re- 
lation enters  the  eddy  viscosity  as  the  proportionality 
coefficient.  Now  in  order  to  substantiate  it  is  necessary 
to  extend  a correlation  -u^ (o) *i/°^ (-r)>a  over  a certain  time 

^ v 

period  commensurate  with  the  relaxation  frequency  of 

smaller  scales  We  shall  study  the  development  of  that  correlation 
from  the  momentum  equation  (77),  which  is  rewritten  as 


after  substituting  for  the  stress  from  (89).  Here 


UJ  «■  C;  fe  y J 


(120) 


is  a relaxation  frequency  of  rank  -+1 . 

Upon  multiplying  (119)  by  u • ^ (t‘ ,-k)  and  averaging, 

* -V 

we  find 
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The  second  term  on  the  right  hand  side  vanishes  on  account  of 
(70).  After  an  integration,  we  have 


A second  integration  gives 


A double  time  integration  involved  in  the  above  and  was  handled 
like  in  (114). 

On  the  right  hand  side  of  (123),  the  last  term  can  be 
computed  from  writing  an  equation  describing  the  evolution  of 
E?,  similar  to  (98).  As  a result,  we  reduce  that  term  to  the 
same  order  and  form  as  the  left  hand  side  of  (123),  but  at  a factor 
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? 

e smaller.  The  details  of  the  calculations  are  similar  to  those 


tr 

s 


leading  to  (116),  and  w’ll  be  omitted.  Thus,  by  neglecting  such 
a term,  we  find 


W.-yT 


(124) 


2 cj 


(?+•) 


z 


(«■) 


025) 


with 


CO 


7 


(«+') 


where  rr0^  and  have  been  defined  by  (63)  and  (65). 

The  formula  for  the  eddy  viscosity  n^°\  obtained  by 
putting  a • o in  (124),  is  similar  to  the  gravitational  exchange 
function  r in  (117),  except  with  a different  controlling  energy. 


VI I IB  Memory  Chain 

We  note  that  the  transport  functions  T , 1°,  T®,  r°  and 
r“ , as  derived  from  the  single  cascade  method,  depend  on  the 
eddy  viscosity  of  first  rank  n'  = n^.  However,  the  approach  to 
equilibrium  of  ^ 1 ■'  entails  a memory  chain,  as  exemplified  by  the 
formula  (125),  and  rewritten  explicitly  in  terms  of  a spectral 
distribution  of  kinetic  energy  F(k),  defined  by  (105), 
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and  the  chain  continues  Here  r.^(x/k^"'  is  an  eddy  viscosity 
of  the  a-th  rank  in  the  x-space.  The  lower  limiting  wave  number 
k^a according  to  the  definition  (50),  is  written  into  the 
argument- 

In  Ref  5 we  have  discussed  a method  of  cutoff  of  the 
memory  chain  by  molecular  damping,  in  which  the  eddy  relaxation 
frequency  will  be  taken  over  by  a frictional  frequency  at  suffi- 
ciently high  rank  of  the  chain.  The  result  helps  in  determining 
the  cutoff  spectrum  near  the  tail  of  the  dissipation  subrange. 
Since  this  is  not  of  special  interest  in  the  present  work,  we  make 
an  inviscid  approximation,  by  writing 


(127) 
0 28) 


with  the  eddy  viscosity  (128)  set  at  a higher  limiting  wave 
number  than  the  eddy  viscosity  (127),  simplifying  (126)  to 

v - r« 


(129) 


IX.  UNIVERSAL  RANGE  OF  SPECTRUM 


The  universal  range  holds  at  sufficiently  large  wave 
numbers,  such  that 


become  independent  of  k Noting  that 

we  can  rewrite  Eqs.  (102),  (103)  and  (104)  in  new  forms  for 
k = °°,  and  subtract  the  new  forms  from  the  original  equations, 
yielding 


- rv  t°  t ^>j°  w f '«  - r + e i-  f 


(131) 


(i32) 

TI  + a jJ--tr  + . 033) 


with  the  notations 

r=  J • 
' - ■ 

<-  *T  ■ 


Q - rjt-j 

v'  = J f fe  -uo  ) 


t - \J  . 

W * u/ 

(134) 
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IX.  UNIVERSAL  RANGE  OF  SPECTRUM 

The  universal  range  holds  at  sufficiently  large  wave 
numbers,  such  that 

^Jdk'ffk)  , f oik1  Cr(k')  , JLjjk'bifoj 

0 •'o  o 

become  independent  of  k Noting  that 

T (k  * °° ) - o , T^(k~o-J  = ° , 

we  can  rewrite  Eqs.  (102),  (103)  and  (104)  in  new  forms  for 
k = °°,  and  subtract  the  new  forms  from  the  original  equations, 
yielding 


- r°t 

t°  -r  ^J°  + f J=  - r + e + ' 

f 

(131) 

t>  \ j°  = -*r+  t 

(132) 

i 

-i 

a 

-f 

— r-  o 0 

/ 

(133) 

with  the  notations 

r = r 

ooj  y 

J - Je( k**) 

/ = J C k-uo  J 

. 

. 

£ - \ J 

w ^ u/ 

. 

034) 
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In  the  present  quasi-stationary  process,  the  coupled 
integral  equations  do  not  have  the  time  t as  an  explicit 
variable,  with  the  understanding  that  the  spectral  distributions 
F,  G and  H may  vary  slowly  with  time  through  the  physical  para- 
meters e,  e and  t . 

; w 
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X INERTIA  AND  DISSIPATION  SUBRANGES 
(•OR  STABLE  AND  UNSTABLE  CONFIGURATIONS 
U a - 1 and  f I ) 


We  refer  the  configurat  on  of  a heavy  liquid  above  a 
lighter  one  a5  unstable  U - i),  and  the  reverse  configuration, 
e.g.  sea  surface,  as  stable  U - -1).  Their  distinction  lies  in 
the  roles  of  the  gravitational  pull:  in  the  unstable  configura- 

tion, the  g»-avi tat'ona 1 puli  becomes  an  energy  source  for  both 
the  surface  elevation  and  the  velocity  which  endows  it,  while  in 
the  stable  case,  the  gravitational  pull  serves  to  produce  the 
kinetic  energy  at  the  expense  of  the  potential  energy,  so  that  the 
two  energies  must  balance  The  gravitational  pull  is  represented 
by  the  exchange  functions  :-i  and  in  Eqs  (131)  and  (132). 

•»  s 

Since  they  operate  only  at  small  wave  number,  they  become  absent  in 
the  inertia  and  dissipation  subranges,  so  that  the  spectral  laws 
will  be  valid  for  both  stable  and  unstable  configurations. 

Under  those  circumstances,  the  equat.ons  (131)  and  (132) 
govern’ng  tne  spectral  di str iDutions  reduce  to 

vi  + y ' ) J -f  . / j /'(  r/K  j = : + r (134)  I I 

^ * 7')C  ' lr  0 35) 

7 * I 

The  first  of  the  system  (134)  and  (135)  becomes  decoupled 
from  the  second,  and  can  be  >olved  independently  The  flow  of 
energy  between  transfer  and  friction  can  be  best  described  by  a 
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differential  form 


■V->)  d 1 

-X  7 f V -r 


— , 2 y f ^ 


obtained  by  approximating 


r-  j 


The  solution  of  (136)  is 


VM= 

/ 


F = -f  ] YY  yi2)'3^'s 


It  follows,  from  v 1 3 7 ) , the  solut  on  or  (135), 


H - t tL 


1 4 2 £ 


The  gene>ai  solutions  ''37)  and  (138)  cover  subranges 
including  the  three  parameters  , and  we  shail  consider  the 
following  three  sub-anges: 

(a)  Inertia  subrange  with  r r » c. t on  i 1 o,  - - o) 

By  writing  o and  = o.we  reduce  (137)  and  (138)  to 


z > - 3 


~(3  <)  k 


H • •=  (f.  > l)  i 


(b)  Dissipation  sub<ange  with  surface  friction  (.  - o) 
When  the  spectrum  (k)  falls  by  friction,  as 


-4?- 


042) 


the  surface  spectrum  is  dissipated  b>  t,  moecular  diffusivity 
a and  become 


' _ (,43) 

(c)  Dissipation  subrange  without  friction  (,  - o) 

By  writing  o,  we  reduce  the  general  solutions  (37) 


and  (138;  to 

F = 

k~7 

J v J K 

(144) 

H - 

Tj/v)  ■ 

(145) 

in  agreement  with  the  He  e^berg  law  or  (1). 

(dj  Ine-tia  Subrange  without  friction 

The  ntrtia  sub-ange  wthout  friction  reduces  further 
Eqs  (135)  and  136J  to  a decoupled  System 

(146) 


f 

* 


(147) 


which  gives  the  solution. 


Vi 
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■,  ■ — wan— 


(148) 


hi 


(149) 


’n  agreement  wt th  the  Ko'mogoroft  law  of  (l).  it  is  to  remarked 
that  tree  medium  and  a movinq  surface  have  two  different  numerical 
coefficients  in  their  inertia  subrange. 
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Xi  bENERAT I ON  OF  TURBULENCE  BY 
uRAYITATIONAl  INSTABILITY  U 1) 


in  v>ew  of  the  unstable  configuration  of  a heavy  liquid 
above  a lighten  one  (i  e z - 1),  a turbulent  interface  may  be 
generated  by  the  gravitat ’ona'  instability  Unlike  the  stable 
sea  surface,  where  the  effects  of  the  gravi tational  pull  balances 
themselves  m the  kinetic  energy  and  the  potential  energy,  the 
present  unstable  conf lgu^ation  relies  upon  the  transfer  function 
tor  transferring  energy  across  a spectrum,  and  upon  the  gravity 
exchange  function  tor  teed  ny  energy  separately  into  each  spectrum. 
Thereto»e,  in  the  production  subrange  where  the  above  mechanism  is 
in  effect,  the  molecular  diss'patons  can  De  neglected,  reducing 
the  equat  on  ot  ene«g>  balance  to 


y(*/k') 


f=r  t- 1 +■  i' 


i- 


\r-  / 


/ k‘  ) 

/ 


(150) 


(151) 


The  equat’on  ot  potential  energy  is  not  useful  here,  because  there 
will  be  no  balance  between  the  kinet.c  energy  and  the  potential 
energy 


The  flow  or  energy  from  the  gravitational  instability  into 
the  wave  coupling  by  inertia  can  be  more  explic-tly  demonstrated 
from  a differential  form  ot  (150)  and  (151), 

p - y _l2_\  / f * o 

(152) 


klC+H 


* 1 


V,  ^ V 


±L 

4 k. 
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The  j'civ itational  puU  n the  unstable  configurat  on 
provide*  a natural  srj'te  of  energy,  t or  maintaining  both  spectra 
in  their  respective  subrange*  of  production,  and  for  keep.ng  them 


from  being  disintegrated  into  d ;s  potion  Under  those  c.ircum 


control  1 -ng  the  eddy 


stances,  the  vort’C’ty  funct'ons  J 


di?s  pations,  can  be  neglected,  reducing  (152;  and  t < 53 ) to 


wher  e 


we  find  tne  so>ut  on 


» j . ^1 


uyg 

i t ^ iii r r 1 1 ni  t ' 


is  a frictiona  transition  wave  number  separating  the  two 


la)  Non-f » »ct lonal , k k , and  T - o 


xi : sea  suRrAcr  turbulence 


CONTROLLED  By  GRAVITATIONAL  PULL 

in  a stab  e configuration  (_  -1)  lue  the  sea  surface,  the 

effect  of  gravity  is  to  pull  down  the  surtace  elevat  on,  and, 
during  this  course  or  act  on,  to  disperse  the  ^urtace  iiquid,  hereby 
to  ra'se  the  velocity  of  dispersal  in  order  to  subscribe  to  the 
above  mechanism,  the  some  g- av  1 tat  >ona  I e*>.hdnqe  function,  which 
represents  tne  pu ! * , shouid  not  only  piay  the  'Oie  of  building  up 
the  kinetic  energy  as  an  acce lera t ion , but  also  of  depletinq  the 
potential  energy  a>  a stab  ng  force,  with  equal  but  opposite 
amount  in  the  system  of  eguat  uni  tor  the  balance  between  the 
kinetic  energy  and  the  potential  energy  Since  the  spectrum  of 
the  kinetic  energy  s in  the  production  subrange  as  a result  of 
the  grav i tat ; onal  pull,  the  friction  w !i  be  neglected  by  writing 
i • o m Eg.  (131).  in  addition, we  omit  the  mo’ecular  dissipation 
functions  m both  eguat’ons,  reducing  the  system  of  equations 
( 131 ) and  ( 133;  to 


(163) 


0 64) 


in  order  to  discern  the  energy  flows  across  each  of  the 
two  spectra  more  conven’entiy , we  rewrite  (i63)  and  ('64;  in 
their  differential  fo'di 
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0 


(165) 


(fk)lH  ,JT 

V <■ + ys  

c.  i ' 1 d k. 


i'r- 


d k 


(i  k) 2 H 


+ ■»: 


OJ 


/ 


w 

f 


jV-‘ 

* 


- 0 


(166) 


describing  the  energy  flows  across  each  of  the  two  spectra  which 
are  in  two  different  regmes  of  development  The  spectrum  of 
kinetic  energy  is  in  Us  early  stage  of  development  in  k-space, 
in  view  of  its  regime  of  production.  Therefore  its  vorticity  is 
not  .yet  formed,  permuting  the  approximation  J o.  On  the 
other  hand,  tne  spectrum  of  potential  energy  is  in  its  later 
stage  of  development  in  k-space,  in  view  of  its  stabilization  by 
gravity  Therefore  the  large  wave  numbers  controlling  the 
potential  spectrum  bring  its  vorticity  to  saturation,  permitting 
the  approximation  Jw  j^,  and  at  the  same  time  rendering  n‘ 
negligible  As  a result,  the  system  of  equations  (165)  and  (166) 
simplify  to 


(167) 


/V  _ 

4-  1=0  . 

dk 

The  system  of  integr o-differential  equations  (167)  and  (168) 
yield  the  solutions 


(168) 


F = j k~J  069) 

H = F i'5 

t 
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(170) 


XIII  SIMILARITY  THEORY 


Unstable  surfaces  refer  to  turbulent  surfaces  of  a heavy 
liquid  above  a light  one,  and  turbulence  is  generated  as  a result 
of  Taylor  instability.  Stable  surfaces  refer  to  the  surface  of 
a heavy  liquid  below  a light  fluid,  eg.  a sea  surface. 

We  shall  summarize,  in  Table  1,  the  results  of  the 
spectral  distributions  obtained  from  the  present  analytical  theory. 
This  gives  us  an  opportunity  of  outlining  their  fundamental 
mechanisms  and  governing  parameters,  and  of  supplying  an  addition- 
al dimensional  theory  of  the  spectral  laws  For  the  sake  of 
abbreviation,  we  shall  omit  the  surface  tens’on  in  the  dimensional 
considerations  The  numerical  coefficients  will  also  not  be  written. 
We  distinguish  the  following  subranges: 

A Production  by  gravitational  acceleration 

This  subrange  exists  for  unstable  surfaces  and  is  absent 
for  stable  surfaces 

The  trictionless  case  is  governed  by  the  parameter  g, 
and  on  dimensional  considerations  can  find  the  spectra 

f H * i~3  (172) 

tor  the  kinetic  energy  and  the  surface  elevation,  respectively 
The  frict'onal  case  has  a second  parameter,  which  is  the 
frictional  wave  number 


(173) 
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Since  the  eievat’on  of  the  surface  is  opposed  by  the  friction 
to  its  first  power,  we  find 


Y 


P =■ 


H = (k  /k)z i 3 . 


(174) 


B Inertia 

The  ipectra'i  laws  'n  the  inertia  subrange  are  independent 
of  the  gravitational  effects,  and  are,  therefore,  common  to  stable 
and  unstaD'e  surfaces 

The  trictionless  laws  are  governed  by  the  parameters  e 
and  , and  agree  with  the  Ko'mogoroff  laws: 


h 


H 


(175) 


Their  dimens  ona  der’vat'ons  are  well  known  and  will  not  be 
repeated  here 

If  the  friction  > s dominant,  the  kinetic  energy  is  trans- 
ferred across  the  spectrum  to  secure  a balance  between  the 
friction  and  the  saturated  vorticity  The  governing  parameter 
is 


n 


Y 


(176) 


and  a dimensional  considerat  on  with  this  parameter  gives  a 
spectrum 


(177) 
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: 


Sece  the  surface  spectrum  should  be  proportional  to  e 


according  to  (135),  a dimensional  analysis  using  the  parameters 


fl  and  t yie'ds 

* c 


hi 


( t 


( 


/fiY)  k 


-/ 


(178) 


C Fnctionless  eddj  dissipation  by  gravitational  pull 
This  mechanism  of  wave  dissipation  controls  a stable 
surface  only.  The  governing  parameters  are  the  potential 
vorticty  Jw  and  the  gravity  g Since  the  build-up  of  the 
kmetic  energy  occurs  at  the  expense  of  the  potential  energy 
with  a vort ‘City  j , the  spectrum  of  the  kinetic  energy  is 


F 


- J 


By  balancing  the  gravitational  accelerdtion  with  the 
nonlinear  eddy  transfer  n Eg  ( ] 67 ) » we  find 


Upon  substitufng  i'79),  we  find,  f'om  (180),  the  surface 
spectrum 


) 


on  a dtmen-  unai  argument 


D Molecular  dissipation 


The  fnctionless  laws 


agree  with  the  Heisenberg  theory  for  viscous  dissipation. 

in  order  to  transform  those  viscous  laws  into  inviscid 


but  frictional,  laws,  a wave  number 


has  to  be  introduced,  to  convert  the  formulas  ( 1 83 ; and  (184) 


We  see  that  we  must  hove  m 4 and  n • 2 tor  the  formulas 


(185)  and  ()84)  to  be  nvscd,  entailing 


We  conclude  that  tne  above  dimensional  considerations 


enable  the  reproduct. on  of  the  results  ot  the  analytical  theory 


XI V COMPARISON  WITH  OBSERVATIONS 


The  present  cascade  theory  has  predicted  a wavenumber 
spectrum  k'5  tor  the  fluctuations  of  the  surface  in  the 
gravitational  subrange  (HO).  This  law  differs  from  the 
dimensional  law 


H(k)  - "88> 

4 „ 

as  proposed  by  Phillips  for  sea-surface. 

By  introducing  a streaming  velocity  u$  of  the  surface 
as  a reference  velocity,  we  can  make  a change  of  variables  from 
k into  a new  variable 


-n  = 


089) 


having  the  dimension  of  a frequency,  and  convert  (170)  into  a 
spectral  distribution 


H(*)  = 


normalized  to 


1 

0 


H in) 


with 


i / 2 

^7 


/ 


(190) 


(191) 


(192) 
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The  analytical  law  or  spectrum  (170)  is  now  brought 
close  to  the  dimensional  formula  However,  this  similarity  is 
only  apparent,  because  (170)  expresses  a wavenumber  spectrum, 
in  spue  of  the  presence  of  n having  the  dimension  of  a fre- 
quency, while  (4)  was  proposed  as  a rrequency  law.  Moreover, 
the  dimensionless  coefficient  6,  which  is  determined  in  the 
analytic  law,  could  not  have  been  determined  m the  dimensional 
law  (4). 

Several  observations  have  oeen  reported4’^'^, 
measurmq  the  spatial  elevation  fluctuations  as  a function 
of  time  as  they  are  conveyed  to  the  point  of  the  probe  by  a 
streaming  velocity  u.  If  u_.  is  large  in  comparison  with  the 
phase  velocity,  c(k  km),  of  the  elevation  fluctuations  in  the 
gravitational  subrange,  their  spatial  scales  past  the  point  of 
measurements  will  be  observed  as  fluctuations  in  time,  with  a 
frequency  (189),  according  to  the  Taylor  hypothesis.  Here  km 
is  the  wavenumber  corresponding  to  the  peak  of  the  spectrum 
preceding  the  qr avi tational  Subrange.  Although  the  average 
phase  velocity  can  be  Sufficiently  large,  mvaiidatina  the 
Taylor  hypothesis,  the  small  scales  of  the  gravitational  sub- 
range may  qive  a small  enough  c(k  k ) as  to  satisfy  the 

ill 

Taylor  hypothesis,  thus  enabling  the  comparison  between  measure- 
ments and  the  theoretical  prediction  oased  open  (190)  and  (192). 

it  will  be  convenient  to  express  J in  (192)  in  terms  of  the 

W 

mean  square  slope  J , which  is  a measurable  quantity.  For  this 
purpose  and  with  the  use  of  (120).  we  rewrite  (118)  in  the  form 


which,  thfoui]'-  (130),  is  transformed  into 


H 


and  is  inteqrated  to  qive 


1 ! oik  k i 

v 

0 


Ik  F 


The  kinetic  enerqy  spectrum  derives  its  enerqy  mostly  from  its 
production  subranqe,  with  a spectrum  (169).  Therefore,  without 
much  error,  we  can  substitute  for  F from  (169)  with  a cutoff 
wavenumber  kQ  at  the  outerscale  2-/ k0 , reducinq  (195)  to 


or  equivalently 


A further  substitution  into  (192)  yields 


with 


( l 
\ » 


/ * 

/ 


O I 


(194) 


(195) 


(196) 


(197) 


(198) 


(199) 
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8-10 

This  value  is  computed  by  taking  the  empirical  estimates  ” : 

. 

Uy)  = ° Vj  ^ u'  =•.  / 5-  «. 

* 

where  u is  the  frictional  velocity,  and  u1  is  the  standard 

ft 

deviation  of  velocity  fluctuations. 

It  is  to  be  noted  that  the  mean  square  slope  has 
been  measured  experimental ly  and  reported  by  Wu^ . In  Fig.  1 
we  have  plotted  the  experimental  values  of  6 based  upon  the 
spectral  data  of  Volkov^,  and  found  them  in  good  agreement 
with  the  theoretical  predictions  (198)  and  (199). 


(200) 
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XV  CONCLUSIONS 


The  spectral  laws  in  (172)  are  valid  for  the  production 
subrange  in  an  unstable  surface,  and  are  not  \/a  1 i d in  a stable 
surface  as  claimed  by  Phillips,  see  ^ 3 ) and  (4). 

_5 

A stable  surface  has  the  k law  \170)  in  the  eddy 
dissipation  or  gravitational  subrange.  The  same  power  (143) 
holds  in  the  molecular  dissipation  subrange  If  the  coefficients 
in  the  two  formulas  do  not  differ  much,  the  two  laws  of  the  same 
power  will  appear  in  a continuous  succession  This  explains  why 
the  5th  power  law  is  so  easily  found  on  sea-surface  turbulence. 

Under  the  circumstances  where  the  Taylor  hypothesis  is 
valid,  the  analytical  law  can  be  brought  to  a form  (190)  which 
has  an  appearance  analogous  to  the  dimensional  law  (4)  proposed 

4 

by  Phillips  . The  coefficient  6,  which  was  proposed  as  a 
universal  constant  in  the  dimensional  theory,  becomes  a function 
of  the  dimensionless  sea-surface  slope,  as  predicted  by 

(198).  The  agreement  between  the  measured  relation  for  6 
and  the  theoretical  prediction  (198)  is  shown  satisfactory  in 
Fig.  1. 
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Fiqure  Leqend 


Fiq.  1 Variation  of  coefficient  8 in  Eq.  (192)  with  mean 
square  slope  0^.  The  experimental  points  are 
obtained  by  usinq  the  spectral  data  reported  by 
Volkov^.  The  solid  line  represents  the  theoretical 
prediction  of  Eqs.  (198)  and  (199). 


